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Abstract 



We have analyzed the results of the latest terrestrial neutrino oscillation 
experiments in the framework of a model with mixing of three massive neu- 
trinos and a neutrino mass hierarchy {mi <^ m2 <^ ^s)- In this model, 
oscillations of the terrestrial neutrinos are characterized by three parame- 
ters, Am^ = mg — mf and the squared moduli of the two mixing matrix 
elements Ues and Using the results of disappearance experiments and 

solar neutrino experiments, it is shown that only two regions of possible values 

2 ii2 ii2 ii2 

of {Uesl I^msI allowed: I. {Uesl l^/^sl both small and II. 

I 1 2 I 1 2 

{Uesl is small and |C/^3| is large. If the mixing parameters are in the region I, 
t'/i ^ i^e oscillations are suppressed. In this case the LSND indication in favor 
of ^ z^e oscillations is not compatible with the negative results of all other 
experiments. If the mixing parameters are in the region II, ^ z^,- oscilla- 
tions are strongly suppressed. If massive neutrinos are Majorana particles, 
our analysis shows that neutrinoless double-beta decay could be observed in 
the experiments of the next generation only if the mixing parameters are in 
the region I. 
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I. INTRODUCTION 



The problem of neutrino masses and mixing is the most important problem in neutrino 
physics today. The search for the effects of neutrino masses and mixing is generally con- 
sidered as one of the important and fertile ways to probe new scales in physics. As is well 
known, neutrino masses and mixing naturally arise in GUT models. The see-saw mechanism 
lU of Majorana neutrino mass generation, that is characteristic of GUT models, is the only 
known mechanism that allows one to explain in a natural way the smallness of neutrino 
masses with respect to the masses of all the other fundamental fermions. 

At present, there are several indications in favor of neutrino masses and mixing. It is 
a general opinion (see, for example, Ref. p|) that the most convincing indications in favor 
of a small neutrino mass difference (Am^ ~ 10~^ eV^) and neutrino mixing come from the 
results of the solar neutrino experiments. In all four solar neutrino experiments (Homestake, 
Kamiokande, GALLEX and SAGE) ||^ the observed event rate is significantly lower than 
the event rate predicted by the Standard Solar Model (see Refs. A phenomenological 

analysis of the data that does not depend on the predictions of the Standard Solar Model 
indicates that the data of the different experiments are not compatible with each other if 
there are no transitions of solar z/g's into other states 0. 

In a recent experiment [^, the LSND collaboration found a positive indication in favor 
of z/g oscillations (see, however, also Ref. 0). The LSND collaboration reported the 

detection of nine UeP — ^ e~^n events at a distance of about 30 m from a target in which 
neutrinos are produced in the decays of stopped tt+'s and yU"*"'s. The LSND collaboration 
calculated p| a background of 2.1 ±0.3 events. As shown in Ref. [§], there is a region in the 
plane of the parameters Am^ and sin^ 29 {6 is the mixing angle) in which the LSND result 
is compatible with the exclusion plots obtained in other experiments [[l0| , [ll| searching for 
transitions. 

In this paper we present an analysis of the results of the latest disappearance and ap- 
pearance terrestrial neutrino oscillation experiments in the framework of a general scheme 
in which flavor neutrino fields are mixtures of three neutrino fields with masses that satisfy 
a hierarchy relation. In this scheme the oscillations among all active neutrinos (z/^ ^ u^, 
'^fj, ^ i^T, ^ are rather strongly constrained. We will show that, in the region of the 
values of the mixing parameters that is allowed by the results of reactor and accelerator 
disappearance experiments, the positive LSND result is not compatible with the negative 
results of all other experiments if there is a hierarchy of couplings in the lepton sector. 

From the results of LEP experiments which measured the invisible width of the Z boson 
(see Ref. [|12|) it follows that only three flavor neutrinos exist in nature. However, LEP 
data do not provide a constraint on the number of massive light neutrinos. The number of 
massive neutrinos depends on the neutrino mixing scheme. If the total lepton number is 
conserved, neutrinos are Dirac particles and the number of massive neutrinos is equal to the 
number of neutrino flavors. If the total lepton number is not conserved and only left-handed 
neutrino fields enter in the total Lagrangian (Majorana mass term), neutrinos with definite 
mass are Majorana particles and the number of massive neutrinos is again equal to the 
number of neutrino flavors. In the most general case of neutrino mixing, left-handed and 
right-handed flavor neutrino fields enter in the mass term and the total lepton number is not 
conserved (Dirac and Majorana mass term). In this case massive neutrinos are Majorana 
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particles and the number of particles with definite mass is twice the number of neutrino 
fiavors. In the case of Dirac and Majorana mass term there is a very attractive mechanism 
which explains the smallness of the masses of light neutrinos, i.e. the see-saw mechanism 
HI]. If one assumes that the total lepton number is violated at the GUT scale Mqut by 
a right-handed Majorana mass term and the parameters that characterize the Dirac mass 
term are of the order of the corresponding lepton or quark masses, then in the spectrum of 
Majorana particles there are three very heavy Majorana particles with masses Mi ~ Mqut 
and three light neutrinos with masses rrik — ''^Fk/^GUT (here k = 1,2,3 and mpk is the 
mass of the up-quark or charged lepton in the k^^ generation). In the see-saw scheme it is 
expected that the fields of fiavor neutrinos are predominantly mixtures of the fields of the 
three light Majorana neutrinos. Thus, from a general theory of neutrino mixing it follows 
that the possibility of three massive neutrinos corresponding to the three existing fiavor 
neutrinos is quite plausible. 

It is to be noted that, if the number of massive neutrinos is more than three, oscillations 
between active and sterile neutrinos (neutrinos that do not take part in the standard weak 
interactions) are possible (see, for example, Ref. p!3|JI^ ). Future solar neutrino experiments 



(SNO [0, Super-Kamiokande [jT6[) could allow to check in a model independent way |T7 



if solar B z/g's transform into sterile states. Some possible scenarios with mixing of more 
than three massive neutrinos were considered in Ref. ||18|| . 

In this paper we will assume that three fiavor neutrino fields (with a = e, fi, r) are 
given by a superposition of three (Dirac or Majorana) fields u^l with mass rrik as 



k=l 



Here U is a unitary mixing matrix. 

The masses of all the known fundamental fermions (leptons and quarks) satisfy hierarchy 
relations. It is natural to assume that neutrino masses also satisfy the following hierarchy 

rrii <^m2 <^ms ■ (1.2) 



If neutrino masses are generated with the see-saw mechanism, the hierarchy relation ( p..2|) 
is a consequence of the hierarchy of the lepton or quark masses. 

We will consider oscillations of terrestrial neutrinos under the assumption that the 
squared mass difference Ar/igi = — mf is small and can be relevant for an explanation of 
the deficit of solar u^s observed by all solar neutrino experiments (say Am^i — 10~^ eV^ 
as suggested by the MSW explanation of the solar neutrino problem |1l9| , p0| ). Thus, for 
all experiments with terrestrial neutrinos we have Am2iL/2p -C 1, where L is the distance 
between the neutrino source and detector, and p is the neutrino momentum. 

Under these assumptions the amplitude of z/q, — z/^ transitions can be written in the 
following form 

A.^., = e->^-^ I Yl ^mK\i +Up^K^^w [ -)^ ) \ , (1-3) 

I ||=oo,G 
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TTli 

with Am^ = rrig — ml and = \ + Tn\ ~ p H -. Then, taking into account the fact 

that, due to the unitarity of the mixing matrix, 

UpkK, = - U^sU:, , (1.4) 

k=l,2 

for the probabihty of — >■ z/^j transitions with f3 ^ a we obtain the following expression 
(see, for example, Ref. 

1 , / Am^L 
- Ay .y„ \ 1 — COS 

2 '^"•'^'3 V 2p 



Pu^^up = 7J ( 1 - COS ) , (1.5) 



where 

A..;., = A.,;.. = 4 |?7„3|' |t/;33|' (1-6) 

is the amphtude of z/q, ^ 1^/3 oscillations. 

The expression for the probability of z/^ to survive can be derived from Eq. ( |1.5| ) and the 
conservation of the total probability as 

P,„^,„ = 1 - E ^-'^-^'^/^ = 1-9 ^-'^■^-^ 1 - cos — , (1.7) 

where the oscillation amplitude By^.^^ is given by 

Pyo.\v^ = Y^uc,;up ■ (1-8) 



Using the unitarity of the mixing matrix, from Eqs. (|1.6| ) and (|1.8| ) we obtain 

S..;.. =4|t/„3|'(l-|t/„3|') . (1.9) 

We wish to emphasize the following important features of terrestrial neutrino oscillations 
in the model under consideration. 

1. All oscillation channels [un ^ z/e, z^„ ^ Ur, ^ ^t) are open and the transition 

9 2 2 

probabilities are determined by three parameters. Am , |?7e3| and I^T^s] (from the 
unitarity of the mixing matrix it follows that |?7r3|^ = 1 — \Ue3f — |?7/x3|^)- 

2. The oscillations in all channels are characterized by the same oscillation length 
-^osc = 47rp/Am^ and the amplitudes of the inclusive transitions Ua —>■ Va and ex- 
clusive transitions ^ '^fi {(^ P) are related by the relation ([1.81). 

3. The probabilities of z/^ vp and z/^ — are equal, i.e. 

Pua^UfS = • (-'-•-'-0) 

This relation is a consequence of the fact that the oscillation probabilities depend only 
on the squared moduli of the elements of the mixing matrix. In a general case of mixing 
of three massive neutrinos, Eq. (|1.10|) for a ^ (3 is satisfied only if CP is conserved in 
the lepton sector. Notice that Eq. ( |1.10D for a = /3 is a consequence of CPT invariance. 



The expressions ( [1.5D and ( [1.7|) have the same form as the standard expressions for the 
transition probabihties in the case of mixing between two massive neutrino fields. In the 
latter case only oscillations between two flavors are possible and the oscillation probabilities 
are characterized by two parameters, Am^ and sin^ 26. The important difference between 
the two-neutrino mixing scheme and the scheme with three-neutrino mixing and a mass 
hierarchy that we consider in the present paper is that the second scheme allows simultaneous 
transitions among all three flavor neutrinos. Some preliminary results of the analysis of the 
terrestrial neutrino oscillation experiments in the framework of this scheme were already 
presented in Ref. |^ and in this paper we wish to elaborate the details and present further 
new results. 



II. ANALYSIS OF THE RESULTS OF NEUTRINO OSCILLATION 

EXPERIMENTS 

We will present here an analysis of the results of the latest oscillation experiments with 
terrestrial neutrinos in the framework of the model with mixing of three massive neutrinos 
and a neutrino mass hierarchy. We will also consider the results of the LSND experiment, 
in which positive indications in favor of ^ Ve oscillations were found [H. The new 
experiments CHORUS and NOMAD searching for transitions are under 

way at CERN. These experiments are sensitive to Am^ in a few eV^ range (which is of 
great interest for the problem of dark matter in the universe) and to small values of the 
oscillation amplitude {Ay^.^y^ > 10^"^). We will discuss possible implications of the results 
of these experiments for ^ z/g oscillations after the projected sensitivity will be reached. 
We will also consider possible implications of the results of neutrino oscillation experiments 
for neutrinoless decay ((/3/5)oiy) in the case of massive Majorana neutrinos. 

A. Reactor and accelerator disappearance experiments 

Let us start with the analysis of the results of reactor and accelerator disappearance 
experiments searching for z/g — > and z/^ transitions. No indication in favor of 

neutrino oscillations was found in these experiments. We will use the exclusion plots obtained 



in the Bugey reactor experiment P3 and in the CDHS and CCFR84 accelerator experiments 



26| , |27[| . At fixed values of Am^, the allowed values of the amplitudes B^^.^^ and B^^.^^^ are 
constrained by 

(« = e,/i). (2.1) 

The values of B^, .„ and 5° .„ can be obtained from the corresponding exclusion curves. 
The parameters \Ua3f (with a = e,fi) are expressed in terms of the amplitudes By^.^^ 



as 



|f/„3|' = -(l±v/l-5w.j ■ (2.2) 

It is obvious that, due to the symmetry of Eq. (|1.9| ) with respect to the interchange \Ua2,\^ ^ 
1 — l^asl^, two values of jf/asl^ correspond to the same value of By^-y^. Thus, from the 
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negative results of reactor and accelerator disappearance experiments, one can infer that 
the parameters |?7q,3|^ at fixed values of Am^ must satisfy one of the following inequalities 



f/asi' < al (2.3) 
f/a3|'>l-a°, (2.4) 



<A[^-xl^-B'] ■ (2-5) 



with 

We will consider values of Am^ in the interval, 10~^ eV^ < Am^ < 10^ eV^, that covers 
the range where positive indications in favor of z/^ ^ Ve oscillations were reported by the 
LSND collaboration IQ. The values of the parameters and aP^ in this interval of Am^, 
that were obtained from the exclusion plots of the Bugey experiment |^ and of the CDHS 
and CCFR84 experiments p6| , p7| , respectively, are presented in Fig.|l]. It can be seen from 
Fig.|l] that in the region of Am? under consideration is always small and a° is small for 
Am2 > 0.5 eVl 

From Eqs. (|2.3|) and ( p.4|) it follows that the parameters If/asl^ are either small or large 



(i.e. close to one). However, due to the unitarity of the mixing matrix the parameters jt/esl^ 

I 1 2 I 1 2 I 1 2 

and |[/^3| must satisfy the inequality \Uez\ + \U^z\ < 1 and consequently they cannot be 
both large. Thus, on the basis of an analysis of the results of only the reactor and accelerator 
disappearance experiments, we reach the conclusion that the values of the parameters |f/e3| 
and If/^sl can lie in one of the following three regions: 

I. The region of small \Uez\^ and |f//^3|^. 

I 1 2 I 1 2 

II. The region of small \Uez\ and large |f/^3| . 
III. The region of large l^/esl^ and small |f/^3|^. 



It can be shown [22] that the region III is excluded by the results of solar neutrino 



experiments. In fact, in the model under consideration the survival probability of the solar 
z/g's is given by 



Py.^u. = (1 - If/eal') + \Ue,\' , (2.6) 

(1 2) 

where Pv^-Jve is the survival probability due to the mixing between the first and the second 
generations. If the parameter jf/esl^ is large, from Eq. (|2.(j| ) in the region of Am? under 
consideration we have Py^^y^ > 0.92 for all values of the neutrino energy. This large lower 
bound of the Uf. survival probability is not compatible with the results of solar neutrino 
experiments 0. 

In the following, we will consider, in detail, the region | and the region |l|. We will 
take into account the results of all the latest appearance experiments searching for neutrino 
oscillations. 
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B. The region of small |?7e3p and |f7^3| 



n I 1 2 I 1 2 

In the region H, the parameters \Ue3\ and \U^s\ are small and bounded as 



\U,,f<al and lU^sl' < , (2.7) 

where the quantities and a° are given by Eq.( |2.5| ) and in the interval of Am^ under 
consideration, take the values as presented in Fig.^. The region | is of great theoretical 
interest. It is now well established that the non-diagonal elements of the CKM mixing 
matrix of quarks, V, are small and satisfy the hierarchy |Vi3| -C IV23I <C \Vi2\- Is this feature 
common to the quark and lepton sectors? A hierarchy of couplings in the lepton sector 
similar to the one in the quark sector can be realized only if the values of the parameters 

I 1 2 I 1 2 ITl 

I Ue3 1 and I Ua3 \ lie in the region |. 

I 1 2 I 1 2 in 

If the values of the parameters | Ues \ and 1 1 are in the region |, we can expect that 
^ oscillations are suppressed with respect to z/^ ^ u^. and z/g ^ z/^ oscillations. In 
fact, in the region ||, the amplitude of z/„ ^ oscillations is given by a product of two small 

I 1 2 I 1 2 

quantities If/gsl and jf/^sl , whereas the amplitudes of z/^ ^ z/^ and z/g ^ oscillations are 
linear in one of these two small quantities (see Eq. (|1.6|) ). 



Let us consider, in detail, ^ oscillations. From the results of reactor and accel- 
erator disappearance experiments, we obtain the following upper bound for the oscillation 
amplitude 

A.,;.. <4a°aJ. (2.8) 

In Fig.^], we have plotted the curve that represents this upper bound, which was obtained 
from the results of the Bugey [^], CDHS |]26| and CCFR84 |2^ experiments (the curve 



passing through the circles). In Fig.|^ we have also plotted the exclusion curves obtained in 
the BNL E776 fTg] (dash-dotted line) and KARMEN |TT| (dash-dot-dotted line) experiments 
searching for z^^ z/g transitions. The region allowed by the results of the LSND experiment 
is shown in Fig.^ as the shadowed region between the two solid lines. Taking into account 
that A^^.^^ < By^.y^ (see Eq.( [L.8|) ), we have also plotted, in Fig.||, the exclusion curve for 
By^-y^ found in the Bugey experiment (dashed line). It can be seen from the figure that for 
small values of Am^ (Am^ < 0.5 eV^) this bound on Ay^.y^ is stronger than the direct bound 
obtained by the BNL E776 and KARMEN experiments. It is also seen from the figure that 
this bound is not compatible with the result of the LSND experiment for Am^ < 0.2 eV^. 

Fig.0 shows that, in the range of Am^ under consideration, with the exception of the 
interval 10 eV^ < Am^ < 60 eV^, the limits on the oscillation amplitude Ay^.y^ that can be 
obtained from the results of disappearance experiments are more stringent than the limits 
obtained from the experiments searching for z/^ — > z/g transitions. From Fig.|^, it can be 
seen that the LSND result is not compatible with the results of the reactor and accelerator 
disappearance experiments in the range of Am^ under consideration, with the exception 
of the interval 5eV^ < Am^ < 70 eV^. In this interval of Am^ the result of the LSND 
experiment is only marginally compatible with the exclusion curve obtained in the BNL 
E776 experiment. 

Thus, from the results of disappearance experiments we have obtained rather strong 
limits on the allowed values of the oscillation amplitude Ay^-y^ in the region | of the param- 
eters \Ue3f and It^^sl^. We will take now into account also the results of the FNAL E531 
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experiment on the search for z/^ v^- and Vt transitions and the recent results 

of the CCFR95 experiment on the search for —* Uj. transitions. From the exclusion 
plots deduced from these two experiments, for a fixed value of Am^, the values of A^, .^^ and 
^i/ei^r constrained by 

< , (2.9) 

A..;.. (2.10) 

From Eq.( |1.6|) it follows that in the linear approximation in the small quantities \Ue3f 



\U,sf^\A,^.,,^, (2.11) 



2 

and |f/^3| we have 



|2 



^e3| ^\a,,,,^. (2.12) 



Thus, the parameters {U^^f and \Ue3f are determined, respectively, by the amplitudes of 
^ and ^ i^r oscillations. With the help of Eq. (|2.11| ), from the results of the FNAL 
E531 and CCFR95 experiments it is possible to obtain, for Am^ > 4eV^, more stringent 
limits on the value of \Uu3\ than those obtained from the results of the CDHS and CCFR84 

I 1 2 

disappearance experiments. Combining these limits with the limits on \Ue3\ obtained from 
the results of the Bugey disappearance experiment, we have found rather stringent limits on 
the value of the amplitude Ay^.^^. In fact, from Eqs.([^.7D and ( |2.11| ) we have 



<ae<;.. • (2.13) 

The bound (|2.13|) on the amplitude A^^.^^ obtained from the results of the Bugey, FNAL E531 
and CCFR95 experiments is presented in Fig.|^ (the curve passing through the triangles). It 
can be seen in Fig.^ that the results of reactor and accelerator disappearance experiments, 
together with the results of —>■ v-r appearance experiments exclude all the region of the 
parameters Am? and Ay^.^^^ that is allowed by the LSND experiment. It is to be emphasized 

1 2 1 2 

that this result has been derived under the assumption that the parameters \Ue3\ and \U^3\ 
are both small. 

In Fig.pl, we have also shown the region in the Ay .y^-Am? plane that could be explored 



when the projected sensitivity of the CHORUS and NOMAD pj] experiments, which 
are searching for i/^ — >■ transitions, is reached (the region delimited by the line passing 
through the squares). 

From Eqs.( f^.lID and ( |2.12| ) we obtain the following inequality 

Ay .y < -Al Al . (2.14) 

The limits obtained, by using this inequality, from the results of the FNAL E531 and 
CCFR95 experiments on the search for z/^ — > transitions and from the results of the 
FNAL E531 experiment on the search for v^. — > transitions are presented in Fig.^ (the 
dotted line). 
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Finally, since in the region |the amplitude A^^-^^ is very small, we obtain, from Eq. (|1.8| ), 
the following relations between the amplitudes of inclusive and exclusive transitions 

i?..;.. ^ , (2.15) 

5..;.. ^ A..;.. . (2.16) 

A test of these relations will allow us to check if there is a hierarchy of masses and couplings 
in the lepton sector. 



C. The region of small |J7e3|^ and large |f/^3|^ 

We will now consider in detail the region |I| in which the parameter \Ue3f is small and 
{Ufj^sf is large (close to one), i.e. 

|f/e3|'<a° and |U^3|'>l-a° (2.17) 



with the quantities and given by Eq. 

In this region, ^ oscillations are strongly suppressed with respect to t'^ ^ and 
^ oscillations. This is due to the fact that the amplitude A,^^.,^^ is quadratic in the 

small quantities jf/esl^ and (1 — I^T^al^), whereas the amplitudes A^^-^^ and A^^-^^ are linear 

(see Eq.(|L^)). 

Taking into account the unitarity bound {Ue^f < 1 — |t^^3|^, it follows from Eq.( p.l7D 
that the parameter \Ue3f must also satisfy the following inequality 

\Ue3f < al . (2.18) 



From Eqs.( p^.l7D and ( |2.18| ), we obtain the following upper bound for the amplitude of 
z/g Ur oscillations 

A,^;,, <4Min[a°,aJ] aj . (2.19) 

From this inequality it follows that z/g ^ u^- oscillations are strongly suppressed in the 
region y. In Fig.^, we have plotted the upper bound for the amplitude A^^.^^ obtained from 
Eq. ( p.l9|) using the results of the Bugey, CDHS and CCFR84 experiments (the curve passing 
through the triangles). From this figure it can be seen that in all the considered range of 
Am^ the upper bound for A^^-^^^ is very small, varying from about 10"'' to about 4 x 10~^. 

1 2 1 2 

Up to now we have used the limits on the parameters \Ue3\ and \U^3\ that were obtained 
from the results of disappearance experiments. We will now take into account also the 
results of the BNL E776 experiment that has searched for — > Ve transitions. From 
the exclusion plot obtained in this experiment, it follows that, for a fixed value of Am^, the 
value of the z/^ ^ Ve oscillation amplitude is constrained by 

A.,;., < Al^.^,^ . (2.20) 
Since in the region O under consideration, the parameter jf/usl^ is close to one in the 

— I i2 / I i2\ 

linear approximation in the small quantities jt/esl and (1 — \U^3\ j, we have 
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^4|f/,3|' . (2.21) 
From this relation and Eg. (|2.20|) we have 

\Ue3f < . (2.22) 

The upper bound for \Ue3\ obtained from the resuhs of the BNL E776 experiment are 
presented in Fig.^ (the dashed curve). From this figure it can be seen that for 0.7 eV^ < 
Am^ < 10^ eV^ the hmits on the parameter \Ue3f that can be found from the results of the 
BNL E776 experiment are much more stringent than those obtained from the exclusion plot 
of the Bugey experiment (dash-dotted curve) and from the results of the CDHS and CCFR 
experiments with Eg. ( p. 181) (dash-dot-dotted curve). 



From the positive results of the LSND experiment, we can find a region of allowed Yahies 
of {Uesf in the lUe^f-Am'^ plane. In fact, for a fixed value of Am^, from the LSND allowed 
region in the Aiy^.^^-Am'^ plane, we have 



Ai-l <A^.^ < Ai+l . (2.23) 
Hence, with Eg. ( p.21| ) we obtain 



^(-) ^(+) 

^ < < . (2.24) 

The corresponding allowed region is presented in Fig.^ (the shadowed region between the 
solid lines). 

From the ineguality ( |2.22| ), it is possible to derive the following upper bound for the 



amplitude of z/g ^ oscillations 

^.e;..<<;.,«r (2.25) 

In Fig.|^ we have plotted the corresponding boundary curve obtained from the results of the 
CDHS, CCFR84 and BNL E776 experiments (the curve passing through the sguares). From 
this figure it can be seen that for Am^ > 1 eV^ the upper bound on the oscillation amplitude 
^ue\ur varies from 2 x 10^^ to 5 x lO^'' and is more stringent than that obtained from the 
results of disappearance experiments using Eg.( |2.19| ). 

From the results of the experiments searching for transitions, we can obtain 

even stronger limits on the amplitude of v^. ^ Uj. oscillations. In fact, in the region |I| 
under consideration, where It^ual^ — 1, the negative results of the experiments searching for 

I 1 2 

i^T transitions give the following upper bound on the value of \Ur3\ 

l^rsl' < . (2.26) 



From Egs.( p.22D and ( p.26| ) we have the following constraint on the amplitude of z/g ^ Ur 
oscillations 
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A < 







4 



(2.27) 



The corresponding boundary curve obtained from the results of the BNL E776, FNAL E531 
and CCFR95 experiments is presented in Fig.|^ (the curve passing through the circles). 
From this figure it can be seen that for Am? > 4 eV^ the amplitude of ^ Uj- oscillations 
is extremely small (less than 3 x 10~^). 

Furthermore, from Eqs.( p.22D and (|2.26| ), we have 



/i3 1 



< 



+ Al 



(2.28) 



The upper bound for 1 — |f/^3| obtained with Eq. (|2.28| ) from the results of the BNL E776, 
FNAL E531 and CCFR95 experiments is presented in Fig.^ (the dashed line). The solid 
line in Fig.^ represents the upper bound ( p.l7|) obtained from the results of the CDHS and 
CCFR84 experiments. This figure shows that the inequality ( p.28|) enables us to obtain, in 
a wide range of values of Am^, more stringent limits on the values of the parameter It/^xsl^ 
from the results of the experiments on the search for v^, and — ^ transitions than 

those obtained from the results of disappearance experiments. 

Finally, the fact that the amplitude Ay^.^y^ is very small leads to the following relations 
between the amplitudes of inclusive z/^ — > z/^ and exclusive ^ {(3 ^ a) transitions 



R — 4 -1-4 



(2.29) 
(2.30) 



Thus, in the region |T| the amplitudes B^^.^^ and B^^.^^ are determined by the amplitudes 
Ay^-y^ and Ay^.^^. From the results of the LSND experiment and using Eqs. (p.9|) , ( p. 201) , 
(^), (g;^ and (P3D| ), we obtain 



^tL < i?..;.. < Min A(+.l,A 



<k < i?..;.. < Min AW^,A°^^,^ 



A 



(2.31) 
(2.32) 



A test of these relations in future reactor and accelerator experiments could allow us to 
check the indications in favor of ^ i^e oscillations found by the LSND collaboration. 



III. NEUTRINOLESS f3f3 DECAY 

As we have shown in Section [II A] , in the framework of the model under consideration, the 
results of the solar neutrino experiments and those of the reactor and accelerator neutrino 
oscillation disappearance experiments indicate that the parameter \Ue^\ is small. If massive 
neutrinos are Majorana particles, this fact can have important consequences for {PP)ou decay 
experiments. The matrix element of {PP)ou decay is proportional to (m) = '^^U^-rrii (see. 



for example, Refs. [|T^,Q). The results of the experiments searching for (/3/3)oi/ decay can be 
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summarized as \ {m)\ < 1 eV (see Ref. [0). The expected sensitivity of the next generation 



of experiments is \{m) \ ~ 10 ""^ eV [31 



In the framework of the model under consideration, the hierarchy of neutrino masses 



implies that |32 



|(m)| ~ lUesf VAm^ . (3.1) 



In Fig.|^ we have plotted the boundary curve for | (m) | obtained with Eq. ( |3.1[ ) from the 
results of the Bugey experiment in the interval 10~^ eV^ < Am^ < 10^ eV^ (the dash-dotted 
line). From this figure it can be seen that for Am^ < 5eV^, the upper bound for is 
less than the projected sensitivity of the {PP)ou decay experiments of the next generation. 

In the region 0, where {Uesf is small and I^T^s]^ is large, the value of \Ue3f is con- 
strained not only by Eq.( |2.17| ), but also by Eq. (|2.18| ), which takes into account the unitarity 
constraint. In Fig.^ we have plotted the corresponding curve obtained from the results of 
the CDHS and CCFR84 experiments (the dash-dot-dotted line). Furthermore, in region 
m the value of \Ue3f is severely constrained by the results of the BNL E776 experiment 
that searched for z/^ — > transitions (see Eq.( p.22|) ). The corresponding boundary curve is 
presented in Fig.|^ (the dashed curve). From this figure it can be seen that \ {m)\< IQ-^eV 
practically in all the considered range of Am . Thus, if the parameters jf/esl and If/^^sl are 
in the region ||, the observation of {(313)0^ decay becomes a formidable, if not impossible, 
task. 

Finally, from the results of the LSND experiment, for each value of Aim?, an allowed 
range of jf/esl^ can be determined from Eg. ( p. 241) in region In Fig.|^ we have also plotted 



the corresponding allowed region in the Aim?-\ (m) \ plane (the shadowed region between the 
two solid lines). 



IV. CONCLUSIONS 

We have considered a scheme with mixing of three massive neutrino fields and a hierarchy 
of neutrino masses. We have assumed that the squared mass difference — ml is small and 
can be relevant for the suppression of solar z/g's. In this scheme the oscillations of terrestrial 
neutrinos in short-baseline experiments are determined by the values of three parameters, 
one squared mass difference Am^ = m^ — m\ and the squared moduli of the two mixing 
matrix elements f/es and U^^,. 

After the calibration of the GALLEX detector with a radioactive source |Q the indica- 



tions in favor of neutrino oscillations coming from solar neutrino experiments have become 
more significant. For this reason, we believe that the model considered here is the simplest 
and most realistic model of neutrino mixing. It seems very appropriate to analyze in the 
framework of this model all the data from the existing experiments that have searched for 
neutrino oscillations and to infer predictions for the results of future experiments (see also 



Refs. |3|,35|). In our discussion, we have also taken into account the positive indications in 
favor of ^ v^, oscillations that were found in the recent LSND experiment. 

We have shown, by using the results of reactor and accelerator disappearance experi- 
ments, that in a wide range of Am^ (0.5 eV^ < Am^ < 10^ eV^) the parameters |?7e3|^ and 
I^^AtsI^ can be either very small or very large (close to 1). From the unitarity of the mixing 
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matrix and from the results of the solar neutrino experiments it follows that only two regions 
for the parameters {Ue^f and If/^^sl^ are allowed: 

I. The region of small \Ue3f and If/^^sl^. 
II. The region of small \Ue3f and large |f/^3|^. 

The overall situation for the allowed regions in the present theoretical framework is 
summarized in Figs.|^ and M. In these figures the various constraints and the allowed regions 
for the values of the parameters jf/esl and |f/^3| are displayed at the representative value 
Am^ = 6eV^ 

If the parameters \Ue3f and IfZ/^sl^ are in the region |, z/^ ^ z/g oscillations are strongly 
suppressed. We have shown that in the region | the indications in favor of ^ z/g oscillations 
found in the LSND experiment are not compatible with the results of all other experiments 
which have not found any evidence of neutrino oscillations. It is to be emphasized that a 
hierarchy of couplings in the lepton sector analogous to that of the quark sector is possible 
only if the parameters | Ue3 \ ^ and 1 11^3 \ ^ lie in region |. 

2 I 1 2 I -• 

If the parameters \Ue3\ and |?7^3| are in the region |I|, the amplitude of Ue ^ i^t oscil- 
lations is very small in the wide range of Am^ under consideration. In this case there is 
no constraint on the amplitudes of z/^ ^ z/g and z/^ ^ z/^ oscillations and the amplitudes 
that characterize the survival probabilities of z/g's and z/^'s are determined by the ampli- 
tudes A^y^-i,^ and A^^-i^^. The indication in favor of z^^ ^ Ue oscillations found in the LSND 
experiment can be checked in future disappearance experiments. 

We have also discussed the implications of our analysis for (/9/9)oj/ decay if massive neu- 

I 2 I 2 

trinos turn out to be Majorana particles. If the parameters If/gsl and jf/^sl lie in the region 
I and Am^ > 5eV^, {PP)ou decay could be observed in the next generation of experiments. 
If the parameters \Ue3f and l^/^sl^ are in the region |I|, the sensitivity of the {P/3)ou decay 
experiments of the next generation is not sufficient to observe this process. 
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FIGURES 



FIG. 1. Values of the parameters and (see Eq.(p.5|)) obtained from the results of reactor 
and accelerator disappearance experiments for Am^ in the range 10"^ eV^ < Am2 < 10^ eV^ 

FIG. 2. Exclusion regions in the A^^-^^-Am? plane for small |?7e3p and iC/^^sp. The regions 
excluded by the BNL E776 and KARMEN ^ appearance experiments are bounded by 
the dash-dotted and dash-dot-dotted curves, respectively. The dashed line represents the results 
of the Bugey experiment. The curve passing through the circles is obtained from the results 
of the Bugey, CDHS and CCFR84 experiments using Eq. (|2.8| ). The curve passing through the 
triangles is obtained from the results of the Bugey, FNAL E531 and CCFR95 experiments using 



Eq.( |2.13D . The dotted line is obtained from the results of the FNAL E531 and CCFR95 experiments 



usmg Eq.(|2l^). The line passing through the squares bounds the region that will be explored by 



CHORUS and NOMAD. The region allowed by the LSND experiment is also shown (the shadowed 
region limited by the two solid curves). 

FIG. 3. Exclusion regions in the A^^-^^-Am? plane for small jUesI^ and large |C/^3|^. The solid 
line represents the results of the FNAL E531 experiment. The dashed line was obtained from the 
results of the Bugey experiment. The curve passing through the triangles is obtained from the 
results of the Bugey, CDHS and CCFR84 experiments using Eq.( p.l9 ). The curve passing through 



the squares is obtained from the results of the CDHS, CCFR84 and BNL E776 experiments using 
Eq.(2.25). The curve passing through the circles is obtained from the results of the BNL E776, 



FNAL E531 and CCFR95 experiments using Eq.(^37|) 



I 1 2 I 1 2 I 1 2 

FIG. 4. Upper bounds for \Ue3\ in the region of small {Uesl and large |f7^3| . The dash-dotted 
and dash-dot-dotted curves are obtained from the results of the Bugey experiment and from the 
results of the CDHS and CCFR84 experiments, respectively (with the help of Eq.(2.1^)). The 



dashed curve is obtained from the results of the BNL E776 experiment using Eq.(|2]2|). The 
shadowed region within the two solid lines is the allowed region obtained from the results of the 
LSND experiment using Eq.( ^.24|) . 

I 1 2 I 1 2 I 1 2 

FIG. 5. Upper bounds for 1— I 1 in the region of small | C/e3 1 and large | 1 . The solid curve 
was obtained from the results of the CDHS and CCFR95 experiments (see Eq.( p.l7|) ). The dashed 
curve was obtained from the results of the BNL E776, FNAL E531 and CCFR95 experiments using 



Eq.dJI). 



FIG. 6. Boundary curve in the Am?-\{m)\ plane for small If^esl^ and |?7^3|^. The curve was 
obtained from the results of the Bugey experiment. 
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FIG. 7. Boundary curves in the Am^-|(m)| plane for small |C/e3p and large |C/^3p. The 
dash-dotted and dash-dot-dotted curves are obtained from the results of the Bugey experiment 
and the results of the CDHS and CCFR84 experiments, respectively. The dashed curve is obtained 
from the results of the BNL E776 experiment. The shadowed region within the two solid lines is 
the allowed region obtained from the results of the LSND experiment. 

FIG. 8. The regions |and|l|of the values of the parameters |C/e3p and iC/^al^ for Am^ = 6eV^. 
The vertical axis has been expanded logarithmically for |C/^3|^ very small and very large (i.e. close 
to one) using the two different coordinates, |C//i3|^ below the horizontal solid line and (1 — |f/;i3p) 
above the horizontal solid line. We have drawn the constraints given by the results of the Bugey 
and CDHS disappearance experiments, which delimit the allowed regions P and p[ The black region 
is excluded by unitarity. 

FIG. 9. Allowed regions for the values of the parameters I C/e3 1 and|[/^3| for Am =6 eV . The 
vertical axis has been expanded logarithmically as in Fig.^. We have drawn the constraints given 
by the results of the Bugey and CDHS disappearance experiments (as in Fig.^) and the constraints 
given by the results of the BNL E776 (f^ — > fg) and FNAL E531 {v^ — > i^r) experiments. The region 
allowed by the LSND experiment is shown as a lightly shadowed band. The darkly shadowed region 
is excluded by unitarity. The region la is the part of region | which is allowed by all experiments, 
except LSND. The region Ha is the part of region || which is allowed by all experiments, except 
LSND, whereas the region lib is the part of region |^ which is allowed by all experiments, including 
LSND. 
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Figure [1] 
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Figure |2| 
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Figure |9| 
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